AC-driven vortices and the Hall effect in a tilted washboard planar pinning potential 
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The Langevin equation for a two-dimensional (2D) nonlinear guided vortex motion in a tilted 
cosine pinning potential in the presence of an ac current is exactly solved in terms of a matrix 
continued fraction at arbitrary value of the Hall effect. The influence of an ac current of arbitrary 
amplitude and frequency on the dc and ac magnetoresistivity tensors is analyzed. The ac current 
density and frequency dependence of the overall shape and the number and position of the Shapiro 
steps on the anisotropic current-voltage characteristics is considered. An influence of a subcritical 
or overcritical dc current on the time-dependent stationary ac longitudinal and transverse resistive 
vortex response (on the frequency of an ac-driving fi) in terms of the nonlinear impedance ten- 
sor Z and a nonlinear ac response at fi-harmonics are studied. New analytical formulas for 2D 
temperature-dependent linear impedance tensor Zl in the presence of a dc current which depend 
on the angle a between the current density vector and the guiding direction of the washboard PPP 
are derived and analyzed. Influence of a-anisotropy and the Hall effect on the nonlinear power 
absorption by vortices is discussed. 

PACS numbers: 74.25.Fy, 74.25. Sv, 74.25. Qt 



I. INTRODUCTION 

It is well known that the mixed-state resistive prop- 
erties of type-II superconductors are determined by the 
dynamics of vortices which in the presence of pinning 
sites may be described as a motion of vortices in some 
pinning potential^. In the simplest case this pinning 
potential is assumed to be periodic in one dimension, 
and temperature-dependent dc current uniaxial pinning 
anisotropy, provoked by such washboard planar pinning 
potential (PPP) recently has been extensively studied 
both theoretically^"^ and experimentally.^"^ Two main 
reasons stimulated these studies. First, in some high- 
Tc superconductors (HTSCs) twins can easily be formed 
during the crystal growth. Second, in layered HTSCs 
the system of intcrlayers between parallel a6-planes can 
be considered as a set of unidirectional planar defects 
which provoke the intrinsic pinning of vortices.™ 

As the pinning force in a PPP is directed perpendicular 
to the washboard channels of the PPP—, the vortices gen- 
erally tend to move along these channels. Such a guided 
motion of vortices in the presence of the Hall effect pro- 
duces anisotropic transport behaviour for which even (-I-) 
and odd (-) (with respect to the magnetic field reversal) 
longitudinal (||) and transverse (_L) dc nonlinear magne- 
toresistivities pjj'^j^ depend substantially on the angle a 
between the dc current density vector j and the direction 
of the PPP channels ("guiding direction"). 

The dc-current nonlinear guiding problem was exactly 
solved recently for the washboard PPP within the frame- 
work of the two-dimensional (2D) single-vortex stochas- 
tic model of anisotropic pinning based on the Fokker- 



Planck equation and rather simple formulas were derived 
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for the dc magnetoresistivities p 

On the other hand, the high-frequency and microwave 
impedance measurements of a mixed state can also give 
information about the flux pinning mechanisms and the 
vortex dynamics. One of the most popular experimental 
methods for the investigation of the vortex dynamics in 
type-II superconductors is the measurement of the com- 
plex ac response in the radiofrequency and microwave 
ranges. When the Lorentz force acting on the vortices 
is alternating, then due to the pinning the ac resistive 
response acquires imaginary (out-of-phase) component. 
Due to this reason measurements of the complex ac re- 
sponse versus frequency uj can give important informa- 
tion on the pinning forces. 

The very early model of Gittleman and RosenbluroiS 
(GR) considered oscillations of damped vortex in a gar- 
monic pinning potential. GR measured the power ab- 
sorption of the vortices in Pbin and NbTa films over 
a wide range of frequencies uj and successfully analyzed 
their data with the simple equation 



r]x 



kpx = Fl, 



(1) 



where x is the vortex displacement, rj is the vortex vis- 
cosity, kp is the pinning constant, and is the Lorentz 
force. From Eq. ([IJ follows that the complex vortex re- 
sistivity pv is 



(2) 



where p/ is the flux- flow resistivity and tUp = kp/rj is the 
depinning frequency. As follows from Eqs. ([T]) and 
pinning forces dominate at low frequencies (w <^ LOp) 
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where is nondissipative, whereas at high frequencies 
(w 3> Wp) frictional forces dominate and the vortex resis- 
tivity is dissipative. 

The experimental success of this very simple model 
stimulated the attempts to use it for the interpretation 
of the data taken in HTCSs, where the effects of ther- 
mal agitation are especially important due to their low 
pinning activation energies and the high temperatures of 
the superconducting state. As the GR model was de- 
veloped for zero temperature and could not account for 
the thermally activated flux flow and creep, which are 
very pronounced in HTCSs, there was a need for a more 
general model for the ac vortex dissipation at different 
temperatures and frequencies. 

In order to fulfill this aim the vortex equation of mo- 
tion ([1]) was supplemented with Langevin force which was 
assumed to be Gaussian white noise with zero mean and 
the cosine periodic pinning potential was used^^~^^ for 
taking into account the possibility of vortex hopping be- 
tween different potential wells. In the limit of small ac 
current (i. e. for a nontilted cosine pinning potential) this 
new equation of motion was approximately solved by a 
continued-fraction expansion^^"^'^ using the analogy be- 
tween a pinned vortex and a Brownian particle motion in 
a periodic potential. As a result, the complex resistivity 
p„ which generalizes the GR's Eq. ([2]) has the form (see 
Eq. (8) in Ref. 11) 

iPv/Pf) = [«(t^/t^o) + i^oo]/[l + i(w/wo)], (3) 

where vqo is a creep factor that grows monotonically with 
temperature increasing from vqq — (no flux creep) to 
'^oo — 1 (flux flow regime) and ojq is a characteristic 
frequency (nonmonotonic in temperature) which, in ab- 
sence of creep, corresponds to the depinning frequency 
ujp. If the frequency u is swept across the temperature- 
dependent frequency wq, the observed Repy increases 
from a low frequency value to the flux- flow value p/ while 
Impi, exhibits a maximum at ujq. Thus, we can sum- 
marize that the temperature-dependent ac-driven vortex 
motion problem has been exactly solved so far only for 
the one-dimensional (ID) nontilted cosine pinning poten- 
tial at a small oscillation amplitude of the vortices. 

At the same time, the examination of a strong ac- 
driving (that is interesting both for theory and for differ- 
ent high-frequency or microwave applications) evidently 
requires to consider strongly tilted pinning potential. Ac- 
tually, if at low temperatures and relatively high frequen- 
cies in nontilted pinning potential each pinned vortex will 
be confined to its pinning potential well during the ac 
period, in the case of strong ac + dc driving current the 
running states of the vortex may appear when it can visit 
several (or many) potential wells during the ac period. 

The aim of this work is to suggest a new theoretical ap- 
proach to the study of temperature-dependent nonlinear 
ac-driven pinning-mediated vortex dynamics based on an 
exact solution (in terms of a matrix continued fraction) of 
the same equation of vortex motion, as was discussed by 
Coffey and Clem (CC) in the seminal paper— (see below 



Eq. Q which has an additional Hall term) . This new ap- 
proach substantially generalizes the CC's results because 
the two-dimensional (2D) Langevin equation for the non- 
linear guided motion in a tilted cosine PPP in the pres- 
ence of a strong ac current at arbitrary value of the Hall 
effect has been exactly solved. For this exact solution 
we used the matrix continued fraction technique earlier 
suggested and later extensively employed for calculation 
of ID nonlinear (ac -|- (ic)-driven response of overdamped 
Josephson junction with noise in Refs. 14, 15. 

As a result, two groups of new findings were obtained. 
First, for previously solved in Refs. 2 and 3 the 2D dc- 
problem of the influence of an ac current on the over- 
all shape and appearance of the Shapiro steps on the 
anisotropic dc pjf'^j^ - CVC's was calculated and ana- 
lyzed. Second, for the ac current at a frequency w plus 
dc bias the 2D nonlinear time-dependent stationary p'^^'^^ 
ac-response on the frequency ui in terms of nonlinear 
impedance tensor Z and a nonlinear ac response at lu- 
harmonics was studied. 

The organization of the paper is as follows. In Sec. II 
we introduce the model and the basic quantities of in- 
terest, namely, the average two-dimensional electric field 
and the Fourier amplitudes for the averaged moments 
(r™). In Sec. HI we present the solution of the recurrence 
equations for the Fourier amplitudes in terms of matrix 
continued fraction and introduce the main anisotropic 
nonlinear component of our theory - the average pinning 
force, divided into three parts. In Sec. IV we discuss the 
w-dependent dc current magnetoresistivity response with 
different (from A to E subsections) aspects of this prob- 
lem. Section V (with subsections from A to H ) repre- 
sents different problems related to nonlinear anisotropic 
stationary ac response. In Sec. VI we conclude with a 
general discussion of our results. 

II. FORMULATION OF THE PROBLEM 

The Langevin equation for a vortex moving with ve- 
locity V in a magnetic field B = nB (i? = |B|, n = nz, 
z is the unit vector in the z direction and n = ±1) has 
the form 

?7V -t-na^v X z = Fi +Fp + Ff,i, (4) 

where F^ = »T.($o/c)j x z is the Lorentz force ($0 is 
the magnetic flux quantum, and c is the speed of light), 
j = j(^) — j''^ + j"''^ cosujt, where i'^'^ and are the dc 
and ac current density amplitudes and w is the angu- 
lar frequency, Fp = —VUp{x) is the anisotropic pinning 
force {Up{x) is the washboard planar pinning potential), 
Fth is the thermal fluctuation force, 77 is the vortex vis- 
cosity, and an is the Hall constant. We assume that the 
fiuctuational force Fth{t) is represented by a Gaussian 
white noise, whose stochastic properties are assigned by 
the relations 

(FthAt)) = 0, {FthAt)Fth^3it')) = 2T?75„<5(t - t'), (5) 
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FIG. 1: System of coordinates xy (with the unit vectors x 
and y) associated with the PPP washboard channels and the 
system of coordinates x'y' associated with the direction of the 
current density vector j; a is the angle between the channels 
of the PPP and j, /3 is the angle between the average velocity 
vector of the vortices (v) and j; Fz, is the Lorentz force; < 
Fp > is the average pinning force, Vlx is the average effective 
motive force for a vortex. Here for simplicity we assume e = 
0. The schematic sample configuration for three cases with 
different values of angle a (the insert): general case, a ^ 
0,7r/2 (a); longitudinal L-geometry, a = 7r/2, j ± y (b); 
transverse T-geometry, a = 0, j _L x (c); in all cases 
and E\\ are transverse and longitudinal (with respect to j- 
direction) electric field components. 



where T is the temperature in energy units, (...) means 
the statistical average, and Fth,i{t) with i — x oi i — y is 
the i component of 'Fth{t)- 

The formal statistical average of Eq. (jj]) is 

ry(v) + naij(v) x z = Fl + (Fp). (6) 

Though (Ff/i) disappears because of the stochastic prop- 
erty in Eq. ([5]), effects of the thermal fluctuation are im- 
plicit in the term (F^) (see below). 

Since the anisotropic pinning potential is assumed to 
depend only on the x coordinate and is assumed to be 
periodic {Up{x) = Up{x + a), where a is the period), 
the pinning force is always directed along the anisotropy 
axis X (with unity vector x, see Fig. 1) so that it has 
no component along the y axis (Fpy = —dUp/dy — 0). 
Thus, Eq. dl]) reduces to the equations 

Vx + Svy = {Flx + Fpx + Fx/)r], 

(7) 

Vy - 6vx = {FLy + Fy)/r], 

where S = ne, e = an /?7, and we omitted index th in 
the Fth for simplicity. Eqs. (I?]) can be rewritten for the 



subsequent analysis in the following form 

Vx = x= {Flx + Fpx + Fx)/ijD, 

(8) 

Vy = y= {pLy +Fy+ SFpx)/iiD. 

where Flx = Flx - SF^y, FLy = F^y + 5Flx, Fx = 
Fx-SFy, Fy = Fy+dFx,D = l + d^ and (F,(t)F, (i')) = 
2TrjDS^.jS{t-t'). 

Our aim now is to obtain from Eqs. ([5]) a rigorous and 
explicit expression of (vx) and (vy) in which effects of the 
pinning and the thermal fluctuation are considered. We 
assume, as usual^'^^"^^, a periodic pinning potential of 
the form Up{x) ~ (C/p/2)(l — coskx)) where k = 27r/a. 
As Fpx = —FpSinkx, where Fp = Upk/2, the first from 
Eqs. dH) has the form 

T{dx/dt) + sin Flx + Fx. (9) 

Here x = kx is the dimensionless vortex coordinate, 
f = rjD/kFp is the relaxation time, F^x — F^x/Fp is 
the dimensionless generalized moving force in the x direc- 
tion. Fx = Fx/Fp, and {Fx{t)Fx{t')) = T6{t - t'), where 
T = 2T/g and g ~ Up/2T is the dimensionless inverse 
temperature. 

Making the transformation x{t) — >• r"^{t) = e~*™^(*) in 
Eq. ^ one obtains a stochastic differential equation with 
a multiplicative noise term, the averaging of which yields 
a system of differential-recurrence relations for the mo- 
ments (r™) = (e^™^) (as described in detail in Ref. [14]), 
viz., 

Td{r"'){t) /dt + [m^/g + impLx] {r"") (t) = 

(10) 

(m/2)((r™-i)(<) - (r™+i)(^)). 

The main quantity of physical interest in our problem 
is the average electric field, induced by the moving vortex 
system, which is given by 

(E) = (n/c)B X (v) = n{B/c){^{vy)^ + («,)y), (11) 

where x and y are the unit vectors in the x and y direc- 
tions, respectively. 

As follows from Eq. ^ 

{vy) =FLy/v + S{vx) (12) 

and so for determination of (E) from Ec^. dTTI) it is suffi- 
cient to calculate the {vx) from Eq. ([9]). This calculation 
gives 

{v.)it) = ^[f^+rcosu;t-{sinx)it)], (13) 
crjL) 

where 

(sinx)(t) = ^[(r)(t)-(r-i)(t)]. (14) 

In Eq. (10) j'^'= = n(jf + 5j^-)/j„ = n{j-- + 6j--)/j,, 
and jc = cFp/^Q. 
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Since we are only concerned with the stationary ac re- 
sponse, which is independent of the initial condition, one 
needs to calculate the solution of Eq. (fTO|) corresponding 
to the stationary case. To accomplish this, one may seek 
all the (r™)(t) in the form 



OO 

k— — oo 



kujt 



(15) 



On substituting Eq. (fTSl) into Eq. (|T0|) we obtain recur- 
rence equations for the Fourier amplitudes FJJ^{uj), i. e.. 



pm+l 



(16) 



where 



Zm,kii^) = 2(j''^ -I- ujTk/m — im/g). 



(17) 



III. THE SOLUTION OF THE PROBLEM IN 
TERMS OF MATRIX CONTINUED FRACTIONS 



The scalar five-term recurrence Eq. can be trans- 
formed into the two uncoupled matrix three-term recur- 
rence relations 



Qm(w)C,„(w) + C„,+i(w) = C.m^i{uj), (m = 1,2,...) 

(18) 

and 



-Q;;(-w)C_„(cj)+C_™+i(cj) = C_™_i(a;), (m = 1,2,, 

(19) 

where Qm is a tridiagonal infinite matrix given by 



j"" z™.o(w) j"^ 

Z,„,i(c^) 

i"" 

; ; -: ■■■) 

(20) 



Q„i(cj) = i 



(the asterisk denotes the complex conjugate) and the in- 
finite column vectors Cm(cj) are defined as 



C„(w) = 



^ 
FV 



and Co 



V ; J 



1 





for 



±1,±2,.. 



(21) 



for m — 0. 



Thus, in order to calculate (sinx)(i) in Eq. ([14]), we 
need to evaluate Ci(aj) and C_i(a;), which contain all 
the Fourier amplitudes of {r){t) and {r~-^){t). Equa- 
tion (|18p can be solved for Ci in terms of matrix contin- 
ued fractions viz., 



Ci(l.) 



-Co, 



(22) 



Qi(^) 



Q2(C^) 



Q3(^) 



where the fraction lines designate the matrix inversions 
and I is the identity matrix of infinite dimension. Having 
determined Ci(a;), it is not necessary to solve Eq. (|19p . 
as all the components of the column vector C_i(a;) can 
be obtained from Eq. (|22p . on noting that 



Foi(c.) = Fo-i*(c.) and F-\uj) ^ FI^Uuj). (23) 

Following the solutions of Eq. and using rela- 

tions (j23p , we can find the dimensionless average pinning 
force {Fpr,){t) (see Eqs. dS])-®, dH]) and ^) which is 
the main anisotropic nonlinear (due to a dependence on 
the ac and dc current input) component of the theory 
under discussion 



•) 



{Fp,){t) = -(sinx)(t) = ^Im(^fee^^-"*), (24) 



fe=0 



where V'o = ^o^('^) ^-nd for fc > 1 we have t/jk 

In fact, Eq. is the expansion of the stationary 

time-dependent (and independent of the initial condi- 
tions) average pinning force {Fpx)(t) into three parts 

{Fj,.,){t) = {Fp^r, -f {Fp,)t, -f (Fp,),^->i. (25) 

In Eq. (gni) (Fpj;)g = -(sinx)^ = ImV^o is the time in- 
dependent (but frequency dependent) static average pin- 
ning force, which will be used for the derivation of the 
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dc magnetoresistivity tensor p^; {Fpx)ti = — (sinx)ti = 
lm('0ie*"*) is the time-dependent dynamic average pin- 
ning force with a frequency w of the ac current input, 
which is responsible for the nonhnear impedance Zi(lj); 
{Fpx)t^^ = -(sinx)t'>^ = Im(?/'fce*''"*) describes a con- 
tribution of the harmonics with k > I into the dynamic 
average pinning force. 



IV. a;-DEPENDENT DC CURRENT 
MAGNETORESISTIVITY RESPONSE 

A. The nonlinear DC resistivity and conductivity 
tensors 

In order to proceed with these calculations we first 
express (see Eq. ([TT]) ') the time independent part of 

{Ey}{t) = {nB/c){vx}{t) as 



v/o 



{nB/c){vx 



{npf/D){j 



dc 



iPf/DHU'y' 



where 



(26) 



l-{^\nxrjf^ = l + {Fpxrjja.. (27) 



In Eq. p6| Pf = B^q/t]c^ is the flux- flow resistiv- 



ity and the Vq can be considered as the (w,j 



dc -ac 



T)- 



dependent effective mobility of the vortex under the in- 
fluence of the dimensionless generalized moving force 
Efx — j'^'^ ill the X direction. In the absence of the ac 
current (see below Eq.(35)) the fg coincides with the 
probability of vortex hopping over the pinning potential 
barrier"^. 

From Eq. (|27p follows another physical interpretation 
of the Vq function, which has a close relationship to the 
average pinning force {Fpx)^ acting on the vortex. Ac- 
tually, it is evident from Eq. (P7| that the {Fpx)^ is con- 
nected to the function in a simple way. 



Then it is easy to show that 



{E.ro = (p//^)[jf(i + '^'(i - ^0 )) - H 



dc'\ 



(28) 



(29) 



From Eqs. ^ and ^ we find the {uj, f T)- 
dependent magnetoresistivity tensor for the dc-measured 
nonlinear law (Eg (w)) = Pgj'^'^ as 



PO ~ \ „dc dc 
\Pyx Pyy 



D 



(30) 



The dc conductivity tensor (Tq , which is the inverse 
tensor to pg , has the form 



1 / 1 



5 



/ dc dc 
[ ^ XX xy 

° " \< O ' PS iD/-o) 



(31) 





ti 



FIG. 2: The dimensionless function VQ{i'^,£,°') numerically 
obtained from Eq. ((27| at g = 30 and Q, = 0.2. 



We see from Eqs. ([30|) and (f3T|) that the off-diagonal 
components of the pg and o-g tensors satisfy the Onsager 
relation {pxy — ^Pyx in the general nonlinear case and 
Cxy = ~'^yx)- All the components of the pg tensor and 
one of the diagonal components of the iTq tensor are func- 
tion of the current density j'*'^, j"'^, and oj through the 
external force F^xi the temperature T, the angle a, and 
the dimensionless Hall parameter 5 = ne. It is impor- 
tant, however, to stress that the off-diagonal components 
of the (Tg are not influenced by a presence of the pinning 
potential barriers. 

In conclusion of this subsection let us consider the lim- 
iting case j"*^ = 0, i. e. we should derive a static current- 
voltage characteristic (CYC). In this case we have from 
Eq. ^ that 



(m + igrWn^ = (.9/2)((r™-^)g - (r™+^)o), (32) 

where the subscript " 0" denotes the statistical average in 
the absence of the ac current. In order to solve Eq. ([5^ 
we introduce, following the calculation of RiskeniS, the 
quantity S„i = {r"^)a/ {r"^~^)o which satisfies next equa- 
tion 



(m + igf')Sra = (5/2)(l - ^™^™+i). 



(33) 



The solution of Eq. p3l) (see details in Ref. [14]) can be 
expressed in terms of the modified Bessel functions I^iz) 
of the first kind of order v (where v may be a complex 
numberJi) as 



(34) 



where p, = igj'''^. Taking into account Eqs. (P7)) . (|34l) . and 
the relation S'l = (r)g = (cosx)g — i(sinx)g we conclude 
that {Fpx)n = ImSi = lra[Ii+p{g)/Ip{g)] and 



1^0 = v^ii"" = 0, w = 0) = 1 + Ini5i/j'*= = 



l + Im[/i+^(.9)//^(g)]/j 



dc 



(35) 
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Note that Eq. (|35|) gives a more simple analytical expres- 
sion for the vq function which was presented in Ref. [2] 
on the basis of a Fokker-Planck approach, namely 

''o'iF) = T^-^ C due~^^Io{2gsm7Tu) (36) 
^ ^ Jo 

with F = where g = Up/2T is the dimension- 

less inverse temperature. In Fig. 2 we plotted i^q {(."^ , S,"") 
graphs at g=30 which demonstrate in the limit of ^'^ = 
the ^'^-dependence for probability of vortex hopping over 
the tilted cosine pinning potential barrier; here and 
are the dimensionless dc and maximal ac current density 
magnitudes (in jc units), respectively (^'^ = j'^/jc,£,°' = 

r/jc). 

B. Longitudinal and transverse DC resistivities 

The experimentally measurable resistive dc responses 
refer to coordinate system tied to the dc current (see 
Fig. 1). The longitudinal and transverse (with respect 
to the dc current direction) components of the electric 
field and Ef, arc related to E^'' = {E^)o and E^'= = 
{Ey)^ by the simple expressions 

^dc ^ ^dc gjj^ ^ _^ ^dc 

(37) 

Ef = -E^" cos a + Ef sin a. 

Then according to Eqs. p7p . the expressions for the 
experimentally observable longitudinal and transverse 
(with respect to the j'*^ direction) magnetoresistivities 
pjjc = £'[j'7j'^ and pf= = Ef/j'^ (where j'^ is the dc cur- 
rent density {j'^Y — iix^Y + ijy'^Y) l^ave the form 
Pf = Pfx sin^ a + pfy cos^ a, 

P± = Pfx sin^ a - P% cos^ a + {p'^'y - pi%) sin a cos a. 

(38) 

Note, however, that the magnitudes of the p'^^'^ and 

p'^f, given by Eqs. (pS)) and applied to the dc current 
responses, in general, depend on the direction of the ex- 
ternal magnetic field B along z axis due to the 5 = ne 
dependence of the function (see Eq. (P7)) ). In order 
to consider only n- independent magnitudes of the p'^^'^ 

and p'^f responses we should introduce the even (-I-) and 
odd (— ) magnetoresistivities with respect to magnetic 
field reversal (p'''^^(n) = [p''-'^{n) ± p''^(— n)]/2) for longi- 
tudinal and transverse dimensional magnetoresistivities, 
which in view of Eqs. (1381) have the form 

pj|^± = (p//i^)[(cos2 a - sin^ a)v'^^ + 

(39) 

+i:>(l± l)sin2 a/2], 
p'^l^ = {pf/D)[Du'^^ smacosa + 5v'^^- 

(40) 

—D{1 ± 1) sinacosa/2]. 



where VQ^{n) = [vo{n) ± t'p (— n)]/2 are the even and 
odd components relative to the magnetic field inversion 
of the function i^Q{n). In the £'[j''j^(j) dependences, which 
follow from Eqs. ([55]) and PO]) . the nonlinear and linear 
(nonzero only for pjl^^ and p"^^) terms separate out in 
a natural way. The physical reason for the appearance 
of linear terms is that in the model under consideration 
for a 7^ there is always the flux-flow regime of vortex 
motion along the channels of the PPP. 

It follows from Eqs. ^ and (gO]) that for a 7^ 0, 7r/2 
the observed resistive response contains not only the or- 
dinary longitudinal p'^^'^'^ (a) and transverse p"^^ (a) mag- 
netoresistivities, but also (as in the absence of ac current, 
see Ref. [3]) two new components, induced by the pinning 
anisotropy: an even transverse p''^~^{a) and the odd lon- 
gitudinal component pj| "^^ («) • 

In the absence of ac current {i/q = I'o (^" = 0,lu = 0)) 
the physical origin of the p''^~^{a) (which is independent 
of e at e ^ 1) is related to the guided vortex motion 
along the channels of the washboard pinning potential in 
the TAFF regime. On the other hand, the p'^^'^^{a) com- 
ponent is proportional to the odd component Vq which 
is zero at e = and has a maximum in the region of the 
nonlinear transition from the TAFF to the FF regime 
at e 7^ (see Figs. 6 and 7 in Ref. [3]). The {f", g) 
dependence of the odd transverse (Hall) resistivity has 
contributions both from the even Uq and from the odd 
Uq components of the voij'^'^, 9) function. Their rela- 
tive magnitudes are determined by the angle a and the 
dimensionless Hall constant e. Note that as the odd lon- 
gitudinal p'^^~ and odd transverse p''^~ magnetoresistiv- 
ities arise by virtue of the Hall effect, their characteristic 
scale is proportional to e ^ 1. 

C. DC response in LT geometries 

In order to analyze the most simple forms of the Py'^j^ 
equations given by formulas (p9)) and pO)) we introduce 
the L and T geometries (see Fig. 1), in which j || x 
(i. e. a = 7r/2) and j _L x (i. e. a and j || y), re- 
spectively. It follows from Eqs. (p9)) and (|40|) that the 
longitudinal py'^" and transverse p^^"*" resistivity for a 
superconductor with uniaxial pinning anisotropy in LT 
geometries vanish (i. e. p'^^^ — p'^'^ = 0) and we obtain 

pt} = <T/D, = ne<j,/i?, (41) 

pt^ = l-eXL/A = ne^o^J^- (42) 

Here <^ = <(j^^ 0., 5), = <(jf , 5), 

If we neglect the Hall terms in Eqs. (|¥T|) and ([^^ . then 
in the absence of an ac current in the L geometry vortex 
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motion takes place along the channels of the washboard 
PPP (the guiding effect), and in the T geometry - trans- 
verse to the washboard channels (the slipping effect) . In 
the L geometry the critical current is equal to zero since 
the FF regime is realized for the guided vortex motion 
along the PPP channels. In the T geometry, i. e. for 
vortex motion transverse to the channels, a pronounced 
nonlinear regime is realized for (7 3> 1 , the onset of which 
corresponds to the crossover point j'^ — jcr, and for 5 ^ 1 
we have jcr = Jo where jc is the critical current. The 
longitudinal even pjj'jt" and transverse odd p'^^Tp resistiv- 
ities are proportional to the even function i'f^rp{^'^, g). In 
the limit j'^, g — > to within terms, proportional to 



^ 1, we have pjj'jt = 1 and p'^j, = ne. The main con- 
tribution to the which is equal to 1 with the same 
accuracy, is due to the guided vortex motion along the 
washboard channels where the pinning is absent. The 
magnitude of p'^'j^ resistivity is described by the Magnus 
force e^'^ which is vanishingly small for a small Hall effect 
for realistically achievable currents j"^ <^ jc/e and the ve- 
locity component (vx) is suppressed, the resistivity p'^J^ 
depends mainly only on the temperature. For g S> 1 the 
/^?L is so small that it cannot be measured (p'J^^ = 
in the limit since e^'^ < 1), and for g ^ 1/2 it 

approaches the value of the Hall constant, e (to within 
terms proportional to <^ 1). 

It is worth noticing that simple Eqs. (|41l) in the T 
geometry allow one to extract from the (^'', ui, g) 
dependences of the measured resistivities pjj^^t and p^^^ 
the dimensionless Hall constant e and the main nonlinear 
component of the model under discussion i^qj^- The latter 
in the absence of ac current, i. e. vq^t, can be used for 
the prediction of the a-dependent pjj'^j'^ resistivities given 
by Eqs. and (gHl) in the case of e < 1. 



„dc- 



D. Guiding of vortices and the Hall effect in 
nonlinear DC+AC regimes 



After derivation of Eqs. (15^ and PO]) let us pro- 
ceed now to a more detailed treatment of the dc vor- 
tex dynamics and the resistive properties associated with 
them in the presence of an ac current. For simplicity we 
will neglect the usually small Hall effect, i. e. we take 
e = 0. As a consequence, the nondiagonal components 
of the dc magnetoresistivity tensor (see Eq. ([501) ) vanish 
{piy = Pyx = 0)- Neglecting the Hall effect, the formu- 
las for the experimentally observed longitudinal p'^^'^ and 

transverse resistivities relative to the dc + ac current 
can be represented as 



pj^ = Pf[{i^S ^ l)sinacosa] — 
= Pfj'^''{Fpx)Q sinacosa. 



(44) 



Therefore, as was pointed out in Ref. [3], even in the 
absence of ac current, under certain conditions in the dc 
current and temperature dependences of the pji^ and p'j^ 
a pronounced nonlinearity appears in the vortex dynam- 
ics and a nonlinear guiding effect may be observed in 
both the inverse temperature g and the current density 
jdc _ ^^d^ consequence of the even parity of Vq 

in and j""" = n^^ (see Eqs. (UHll-dMl) and (gT])) the 
magnetoresistivities and p"^ are even in the magnetic 
field reversal, as they should be neglecting the Hall effect. 

As was shown in Ref. [3], the specifics of anisotropic 
pinning consist in the noncoincidence of the directions of 
the external motive Lorentz force acting on the vor- 
tex, and its velocity (v) (for isotropic pinning || (v) if 
we neglect the Hall effect). The anisotropy of the pinning 
viscosity (which can be defined as the inverse vortex mo- 
bility [{1^)0]^^) along and transverse to the PPP channels 
leads to the result that for those values of {j''''^,g,a) for 
which the component of the vortex velocity perpendicu- 
lar to the PPP channels, {vx)q, is suppressed, a tendency 
appears toward a substantial prevalence of guided vortex 
motion along PPP channels (the guiding effect) over mo- 
tion transverse to the channels (the slipping effect). In 
the experiment, the function 



cot /3 = - ^ = 

pf tana + iy^{^d,g,^^)cota 



(45) 



is used to describe the guiding effect, where (3 is the angle 
between the average vortex velocity vector (v) and the 
current density vector j''^ (see Fig. 1). The guiding effect 
is more stronger when the difference in directions of F^ 
and (v) is larger, i. e., the smaller is the angle /3. Let 
us consider the current and temperature dependence of 
cot /3{^y,g,£_y) for fixed values of the angle a ^ 0,7r/2. 
In the temperature region corresponding to the TAFF 
regime, we have /3 « a and, consequently, at low cur- 
rents guiding arises. At large currents {^y » 1), where 
for vortex motion transverse to the PPP channels the 
FF regime is set up, i. e. the vortex dynamics becomes 
isotropic and we have (v) 1 1 for arbitrary value of the 
angle a. 

Let us now analyze the dc magnetoresistivity depen- 
dences py''^ and p'^j^^, given by Eqs. (HH) and (gni), with 
allowance for the small Hall effect. In this case, the ex- 
pressions for p||'^^, out to terms of order ^ 1, have the 
form 



sm a), 



(43) 



pj|'=+ = p/(i'o + cos^ a + sin^ a), 
p'^+ — p/(i/Q^ — 1) sin a cos a. 



(46) 
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(47) 

p'^ — Pf (<5i^Q ^ + i^Q sin a cos a) , 
Here v^^ are obtained from relations 

i^o" = 1 + (^px)o //" = 1 + ImGo, (48) 
where Go = ^o{j'^'',j'"',^,g)/3'^'' and 

iy:^+ = 1 + ImGj, ly;^- = ImGa . (49) 

In the limit of a small Hall effect (e <C 1) the expres- 
sions for even and odd components of Vq (in terms of G^) 
in the linear approximation in the parameter e tan a ^ 1 
are equal respectively to 

Gj = Go«;J,<^) +ni?+Jtana 

(50) 

Gq = {nRd — Gq)5 ta.na, 

R, EE [d^o/d^i + UVf)id^o/dey)l (51) 

where i/'o = V'o(«Cy j )i J'* and j'' are dc and ac current 
density values, and , are even and odd parts of the 
Rd, respectively. 

As follows from Eq. (|46l) and (|47l) , the behavior of the 
dc current and temperature dependence of pjj'^j^ is com- 
pletely determined by the (j'', glj"^, aj)-behavior of the 
j/q^ dependences. If {.]""/ j"^) ^ 1, i. e. the influence 
of the ac current on the dc response can be considered 
as a small, the linear limit for (j'^'^) dependences, 
following from Eqs. (|46l) and (H7|) . is realized in that 
region of j'^ and g, where (/q = const and = 0, 
while the region of nonlinearity of p'l'^'^^{j'^,g\j°-,uj) de- 
pendences corresponds to those j'^ and g intervals, where 
the dependences Vo'^i^ylg) and i^o^igl^y) are nonlinear. 
Note, that the nonlinearity in the temperature depen- 
dences P||'j'^(5lCy ) can be observed not only at small cur- 
rents (in the TAFF regime), but even at large currents 
> 1 in the case when < 1, where this latter re- 
lation depends on the magnitude of the angle a (for 
cos a < l/^"^ we have < 1 and iy^{g > 1) = 0). 
Thus, the linearity or nonlinearity of the dependences 
pj|°^((7) at dc currents larger than unity depends on the 
magnitude of the angle a. 

E. Shapiro steps and adiabatic DC response 

Before a discussion about the influence of the ac cur- 
rent on the current-voltage characteristic (CVC) of the 
model under discussion it is instructive to consider first 
a simple physical picture of the vortex motion in a tilted 
(due a presence of the dimensionless dc driving force 



< < oo) washboard planar pinning potential (PPP) 
under the influence of the effective dimensionless driving 
force / = Fpx + Flx = - sinx + 

If the temperature is zero, the vortex is at rest with 
= at the bottom of the potential well of the PPP. 
When the PPP is gradually lowered by increasing 
then for < < 1 appears an asymmetry of the left- 
side and right-side potential barriers for a given poten- 
tial well, and in this range of an effective force / 
changes its sign periodically. With gradual ^''-increasing 
there will come a point where ^'^ = 1, and for > 1 
the more lower right-side potential barrier disappears, 
the effective motive force / becomes everywhere along x 
positive and the vortex is in the "running" state, peri- 
odically changing its velocity with a dimensionless fre- 
quency uji = y/ (^'^)^ — 1. So the static CVC of this pe- 
riodic motion at > 1 is a result of time-averaging of 
the stationary time-dependent solution of the equation of 
motion dx/dr = f with t = t/f. Eventually, the proba- 
bility of the vortex overcoming the barriers of the PPP 
i/Q = (^^ = 0,iuj = 0) at zero temperature is 

r 0, ^'^<i, 

1^0 = < (52) 

i. e. the i'o{^'^ > 1) monotonically tends to unity with 
^''-increasing. 

If the temperature is nonzero, a diffusion-like mode ap- 
pears in the vortex motion. At low temperatures (.g 3> 1) 
and < < 1 the thermoactivated flux-flow (TAFF) 
regime of the vortex motion occurs by means of the vor- 
tex hopping between neighboring potential wells of the 
PPP. The intensity of these hops at low temperatures is 
proportional to the exp [— g(l — ^'^)], i. e. strongly in- 
creases with T-increasing and ^''-increasing due to the 
lowering of the right-side potential barriers at their tilt- 
ing. On the other hand, at f '' just above the unity (when 
the running mode is yet weak), the diffusion- like mode 
can strongly increase the average vortex velocity even at 
relatively low temperature due to a strong enhancement 
of the effective diffusion coefficient of an overdamped 
Brownian particle in a tilted PPP near the critical tillji^ 
at ^'' = 1 (see below subsection V. H.). 

Now we consider the influence of a small (^° <C 1) ac 
current density with a frequency uj on the CVC in the 
limit of very small temperatures (g 1). In this case 
the physics of the dc response is quite different depending 
on the ^'' value with respect to the unity. If < 1, the 
vortex mainly (excluding very rare hops to the neighbor- 
ing wells) localized at the bottom of the potential well 
where it experiences a small oj-oscillations. The averag- 
ing of the vortex motion over the period of oscillations 
in this case cannot change the CVC which existed in the 
absence of the ac-drive. 

If, however, ^'' > 1, then the vortex is in a run- 
ning state with the internal frequency of oscillations 
= V (C"*)^ ~ 1- If 7^ the CVC is changed only 
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FIG. 3: The longitudinal CVC C), ^ = 0.2, g = 100, 

showing the Shapiro steps. 



in the second-oder perturbation approach in terms of a 
small parameter ^° <C 1 (as it was shown for the anal- 
ogous resistively shunted Josephson junction problemiS) 
because the CVC is not changed in the linear approxima- 
tion in this case. However, for uj = cji appears a problem 
of a synchronization of the running vortex oscillations at 
the Wj-frequency with the external driving frequency w. 
As a result, the average (over period of oscillations) vor- 
tex velocity is locked in with the w in some interval of the 
dc current density even within the frame of the first- 
order perturbation calculation. The width of this first 
synchronization step (or the so-called "Shapiro step" in 
the resistively shunted Josephson junction problem) has 
been found in Ref. [^°] and the calculation in the spirit 
of this reference gives the boundaries of the where the 
step occurs as 

(d"r/2a<C'<(d+C72C'). (53) 

Here is the current density which gives Ui = 
^yiS.iy - 1 = w, i. e. = Vl +^2. Then the size of 
the first Shapiro step on the CVC is In higher ap- 

proximations (in terms of (C)™)i where m runs through 
all of the integers) the Shapiro steps on the CVC appear 
at the frequencies = muj and fi^ = mw^. The width of 
the m-th step at — ?► is proportional to (^°)™, i. e. 
strongly decreases with m increasing^^. 

In Fig. 3 we plot the longitudinal CVC E'^^';+ {i"^ , i") 
showing the Shapiro steps. The plot in Fig. 3 looks like 
the similar curves discussed earlier^^ for the CVC of the 
microwave driven resistively shunted Josephson junction 
model at r = where the overall shape of the CVC 
and different behaviour of the two types of the Shapiro 
steps in adiabatic limit was explained. Our graph, in 
comparison with the curves of Ref.^^^ is smoothed due to 
the influence of a finite temperature. The longitudinal 
CVC £JJj''ij!'(^'', ^°)-dependence demonstrate several main 
features. First, in the presence of the microwave current 
the dc critical current l^") is a decreasing function of 



the ac driving. The physical reason for such behaviour 
lies in the replacement of the dc critical current by the 
total dc + ac critical current. Second, with gradual 
increasing the zero-voltage step reduces to zero and all 
other steps appear. Such steps are common because they 
do not oscillate and spread over a dc-current range about 
twice the critical current 2^°. These steps are the steps 
of the first kind and they distort the CVC as like as relief 
bump with a concave shift from the ohmic line. With 
further ^"-increasing this relief bump shifts toward higher 
^"^-values. Below this range the steps of the second kind 
appear. These microwave current-induced steps oscillate 
rapidly and stay closely along the ohmic line^S, over a 
dc-current range < ^° — 1. 

To summarize, we can determine three (^'', ^'')-ranges 
where the CVC-behaviour is qualitatively different. For 
large dc bias current densities + 1 < the CVC 
asymptotically approaches the ohmic line without mi- 
crowave induced steps. For an intermediate dc current 
range — 1 < < + 1 CVC curve deviates from 
the ohmic line as a concave bump with the stable steps. 
For lower dc current range < — 1 the steps oscillate 
with microwave current along the ohmic line. With grad- 
ual fi-increasing the size of the steps increases whereas 
their number decreases. 



V. NONLINEAR STATIONARY AC RESPONSE 

A. Derivation of the impedance tensor 

Using Eq. (jlip we determine nonlinear (in the ampli- 
tudes j""^, j'^'^ and the frequency w) stationary ac response 
as 

(E), ^ ((E) (t) - (E)^^) = {nB/c)[{v^)ty-{vy)t^l (54) 

where (E)^ = (nS/c)[-(wj,)^x + {v^)oy] is time- 
independent part of (E)(i) (see also Eqs. ([26|) and ([29]) ). 
whereas {vy)t and {vx)t are time-dependent periodic 
parts of {vy)(t) and {vx){t) which to become zero after 
averaging over a period 2tt/uj of the ac cycle. 
From Eqs. p3|) and ([54| we have 

oo 

{Ey)t - (npfjc/D) ^(j-)'=Rc{Zfe(c.)e*'=-*}, (55) 

k=l 

where 

^fe(w) = <^M-^^fc(w)/(J'^")^ 

(56) 

and 6i^k is Kronecker's delta. 

The dimensionless transformation coefficients Zk in 
Eq. ((55l) have a physical meaning of the A;-th harmonic 
with frequency Clk = kuj in the ac nonlinear {Ey)t re- 
sponse. Equation ([56]) for fc = 1 yields 

^1 = 1 - i^i/f', (57) 
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and using Eqs. (l54l) and (|55p we can express the nonlinear 
stationary ac response on the w-frequency , in terms 
of the nonhnear impedance Zi as 



Ell = {Pf/D){fy^ + 5jr)Re{Zie-*}. 



(58) 



If we put Zi = pi — where pi and Ci are the 
dynamic resistivity and the reactivity, respectively, then 
Eq. ([58|) acquires form 



rpac 



I = iPf/DWy' + 5jT)^pI + CI cosM - ^i), (59) 



where ypf + Ci = l^i| a-^d (/?! = arctan(Ci/pi) are the 
dimensionless amplitude and phase of the ac response on 
the jl% = (j-= + 5jT) cosiot input. 

Similarly, using Eq. (fT2|) and (|54p . we can show that 



= p/jrcoswi- (60) 
and obtain w-frequency ac response as 

= {PflD)Ke{^-'[{D - d'Z.Ux - SZ.jy)]}. (61) 



From Eqs. (|58p and ([61)1 follows that the complex am- 
plitudes of the electric field Ei and the current density 
J = je*"* are connected by the relation Ei = ZJ, where 
Z is the frequency and dc and ac current amplitudes de- 
pendent impedance tensor 



Z{lu) = 





^xy \ 











(62) 



It is relevant to remark the similarity of Eq. (|30l) and 
Eq. (|62p from which follows that for the ac response Zi 
plays the same role as Vq for the dc response. 

However, the connection between Zi and the dynam- 
ical average pinning force {Fpx)ti is more complex than 
the relation between i/^ and {Fp^)'^ (see Eq. (gS])). Tak- 
ing into account the time dependence of the {Fpx)ti, it is 
easy to show that 



{Fpx)ti = Re{^'^'[r{Zi~m. 



(63) 



Equation (j63| gives physical interpretation of the Zi 
impedance and its structure may be compared with 
Eq. (EHD. 

The real quantities Ei = ReEi and j'^'^ — Re J are con- 
nected by the relation Ei = /0°^j, where the ac-response 
resistivity tensor is 



/ o o 

I rxx rxy 

- I ac ac 

\Fyx Hyy 



Note that from Eqs. and (HD) follows that 



(64) 



B. Longitudinal and transverse impedance 
responses 

The experimentally measurable ac resistive responses 
refer to coordinate system to the ac current which di- 
rected, for simplicity, at the same angle a with respect to 
the y axis as the dc current (see Fig. 1). The longitudinal 
and transverse (with respect to the ac current direction) 
components of the electric field E^'^ and E'J^, are related 
to E'^'^{t) and Ey'^it) by the same relations as for the dc 
current (see Eqs. (|57))). The latter is true for the rela- 
tions between the experimentally observable longitudinal 
and transverse (with respect to the j""^ direction) magne- 
toresistivities p'^^'^ = E'^^ j j"- and p°f^ = E°f; j j"- (where 
is amplitude {{ff = {j^")^ + {j^^f). As a resuh we 
have 



= Re{Z||e'"*}, pI^ = Re{Z^e''^*}, 

where 

Z|| = {pf/D)[{D ~6^Zi) sin^ a + Zi cos^ a], 
Z± = {pf/D)[SZi - D{1 - Zi) cos a sin a]. 



(66) 



(67) 



Note, however, that the magnitudes of the and p'^, 
given by Eqs. (|66p. in general (as in the case of dc cur- 
rent), depend on the direction of the external magnetic 
field B along z axis due to the 6 — ne dependence of the 
Zi through the implicit dependence of ipi^Lj) on the j'^'^ 
and j'^'^. In order to consider only n- independent magni- 
tudes of the pif^ and p°j^ resistivities we should introduce 
the even (-1-) and odd (-) longitudinal and transverse 
magnetoresistivities with respect to magnetic field rever- 
sal in the form p'l^''^{n) = [p||^^(n) ± p|f;=^(-n)]/2. 

Let us first separate Zi = 1 — itpi on the even 
Zi{n) = Z^{—n) and the odd Z^ {n) — —Z^{—n) parts. 
If we assume ^pi{n) — il>^{n) + ip^{n), where ^l)^{n) are 
the even and odd parts of il^i{n) (i. e. ''pfin) — [4>i{n) ± 
i/ji{—n)]/2), then we have 

Z+{n) = 1 - inUjy^j^in) - Sjxi^+ (n)]/ {fy - 5^jl), 

Zi{n) = -injc[jyilj+{n) - 5 j xi^i {n)] / [jl - 5'^jl). 

(68) 

From now on, we can present Eqs. (j67|) in the form, 
similar to Eqs. ((39|) and (l40t . with the only difference in 



the change of for Z^ and p'^^f_ for Z^ ^. However 



hereafter it will be suitable for us to present Eqs. (|67| in 
another equivalent form 

Z+ = {pf/D)[{D - S^Z+) sin^ a + Z+ cos^ a], (69) 



{pf/D)Z^ (cos^a- (5^sin^a), 



(70) 



p'"' = Re{i-(a;)e*'^*}. 



(65) 



Zl = {pf/D)[5Z^ 



15(1 - Z+) sin a cos a], (71) 
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Zl = {p}/D){5ZJ + DZ^ sin a cos a) . (72) 



C. The Hall effect and the guiding of vortices in 
nonlinear AC response 

Let us consider peculiarities of the ac resistive re- 
sponses in the investigated model due to the Hall effect. 
Experimentally, three types of measurements of the ob- 
served ac resistive characteristics are possible in a pre- 
scribed geometry defined by a fixed value of the angle a. 
First is ac response measurements which investigate the 
dependence of observed g) resistivities on the 

current density at fixed dc current density Sf^ and tem- 
perature g. Second is the dependence of p|'^^(ff|^", on 

the temperature at fixed and Third, is the depen- 
dence of p'^^^{i'^\i°' ,g) on the dc current density at fixed 
f° and g. The form of these dependences is governed 
by a geometrical factor - the angle a between the direc- 
tions of the current density vector j {t) and the channels 
of the washboard pinning potential. There are two differ- 
ent forms of the dependence of p"^^^ on the angle a (see 
formulas ^^-^^). The first of these is the "tensor" de- 
pendence, also present in the linear regimes (similar to 
the TAFF and FF regimes for formulas and PO]) ). 

which is external to the impedance Zi (see Eqs. (|67l) ). 
The second is through the dependence of Zi on its argu- 
ments £,°'{a) and £f'{a), which in the region of the transi- 
tion from linear in j"^ and j''''^ regimes (at ^"■''^ ^ 1 and 
^o-A ^ 2.) is substantially nonlinear. 

First recall that in the absence of the Hall effect (e = 0) 
there exist only even (with respect to magnetic field in- 
version) impedances j_ - the odd impedances ^_ are 
zero (see Eqs. (l69l) - (l72t V The presence of nonzero value 
of e leads not only to the appearance of a Hall contri- 
bution to the observed ac responses on account of the 
even component Zf of the impedance Zi, but also to 
the appearance of the odd component Z^ , which has a 
maximum in the region of the nonlinear transition from 
the one linear regime (at low and g ^ 1) to another 

linear regime (at large C",^'' and arbitrary g) and is es- 
sentially equal to zero outside this transitional region (see 
Figs. 10, 11 in Ref. [3]). As a consequence, "crossover" 
effects arise: contributions from to effects due to Z^ , 
and vice versa; contributions from Z^ to effects due to 
Z^ . Thus, in the even impedance Z^ (see Eq. (I7T]) ). in 
addition to the main contribution created by the guiding 
of vortices and described by Zf there is present a Hall 
contribution arising due to Zj~. The expression for the 
odd impedance Zj (see Eq. UW ) contains, in addition 
to the Hall term arising due to Z^ , term due to Z^ . 



D. AC response in LT geometries 

In order to study a more simple form of Eqs. ((69)) - 
(17^ we consider first the L and T geometries of the ac 
response (see Fig. 1, insert). 

In L geometry a = tt/2 and jf^ — e^" does not depend 
on n as well as j'j^ = e^"^ for the dc current. As a result, 
-01, L as well as Zi^l = 1— ^V'i.l/Jl'^ have V'j'l — Z^l — 0- 
Finally, from Eqs. follows 

Z+^ = (p//i?)[l + <52(l-Z+J], Z|^^ = 0, (73) 



Z+^ = 0, Z-^ = <5(p;/i?)Z+^. (74) 

If we define p°^^ and C,'^^ as the resistivity and reactivity 
of Zl^j^ impedance, respectively, by the relation Z|^^ = 
~ ^^Ti show that 

p-+ = (p//i?)(l-dm0We), 

(75) 

= -(P//^)eReV'i,L/e. 

Note that experimentally measured quantities l-^i^^l 
and tan(p||^i = I P^l ^-How to obtain and Cfj""/ 
and to compare them with the theoretical formulas (j75p . 
Similar calculations for ZJ ^ = p"'^^ — «C±'^£ yield 

pl^- = n(p;/z?)(6e- Wi,L)/e, 

(76) 

crz = n(p//^)ReV'i,L/e- 

In T geometry (see insert in Fig. 1) a = 0, 2t^ — n^°' 
and j'^ = nS^'^. In this case it follows from Eqs. (P0]) - (P5)) 
that "01, t('^) = —'^i,T{—n), i. e. Vi,t is an odd function 
of n and V'i^t ~ ^- As a result, 

Z+j, = 1 - mV'i^j./C, = 0- (77) 

Then from Eqs. ^ we have 

Z+^^{pf/D)Zl^, Zl^^Q, (78) 



Z-^ = (p;/D)(5Z+j,, Z+^ = 0. (79) 

If Z+y = p'^^+ - and ZJ J, = pj';^, - ^CJ^^;, then 

from Eqs. (ITS]), (HH) we obtain Z J = SZ^rj. and 

pTt^^pI^^ ^T^^^^hT- (80) 

From Eqs. ([50)1 follows that that experimentally mea- 
sured quantities satisfy the simple relations 

I^T tI ~ ^l^lfrl' tan (/sj J, = tan (^ii^j,. (81) 
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At last, from Eqs. dZZ]) and ^ follows that p|['=+ = 
(p//D)p-+ and = (p//^)Ct+, where 

+ = 1 + nIm7/>-T/e, 

(82) 

= nReV-i^x/r, Pi,T = Ci,T = 0- 

E. The power absorption in AC response 

In order to calculate the power absorbed per unit vol- 
ume V (and averaged over the period of an ac cycle) we 
use the standard relation V — (l/2)Re(Ei • J) where 
El and J are the complex amplitudes of the ac electric 
field and current density, respectively. Using Eqs. ([62|) 
and (|67p we can show that 

r = {f/2)p={f/2)ReZ^\. (83) 

After some algebra we obtain that 

{pf/D)[Dsm^a + {l- Dsm^a)RcZi. (84) 

Taking into account that Z\ = \ — iGi, where Gi = 
we conclude that ReZi = 1 — Re(iGi) = 1 + 
ImGi. Then from Eq. (IMl) we have 

p= (p^/i:»)[l + (1 - i:)sin^a)ImGi]. (85) 

In the limit of e ^ 1 we obtain for p a more simple result 

p = p^(l+IniGiCos2a), (86) 

which will be analyzed in detail in Sec. IV. F. 

From Eq. ([5^ follows two simple results for p in LT 
geometries 

PL = {Pf/D){D - S^ReZ^L) = {Pf/D){l - elmV^ii/D, 

(87) 

PT = (p//i^)ReZiT = {pf/D)il + nlmij-j.{n)/n- (88) 

Note that pL and pT in Eqs. ([57)) and ([55]) are equal to 
the expressions for the py'^^ and p||'jt given by Eqs. (|75)) 
and ([82|l . respectively. 

F. AC impedance and power absorption at e <C 1 

Here we analyze the ac + do impedance dependences 
Z^ and Z^ (see Eqs. ^-^), with allowance for the 
small Hall effect (e < 1). In this case Eqs. ([69|) -([72 |) 
become more simple and the expressions for Z^j^ , out of 

terms of order e'^ ^ 1 , have the form 

Z+ ^ pf[Z^ cos^ a + sin^ a), 

(89) 

Z^ — Pf{Z^ — 1) sin a cos a, 



= pfZ cos^ a, 

(90) 

Z^ = pf{SZ^ + Z^ sin a COS a). 
Here Z^ can be obtained from relations 

Zi = Z+ + = 1 - iGx = 

(91) 

1 - i{G+ + Gf) = (1 - iG+) - iG^, 

where Gi = il^iir" , j'^") / f and Z+ = (1 - iGl),Z^ ^ 
—iGi- In the case of a small Hall effect (e <C 1) 
the expression for even (+) and odd (-) components of 
Zi(^''^, (in terms of Gf ) in the linear approximation 
in the parameter e tan a <C 1 are equal, respectively, to 

G+ = Gi«;j,<;|),+ni?+(5tana 

(92) 

Gi (-Gf + ni?^)(5tana, 



i?, ^ [(aV^i/9e,") + {f/n{d^i/d^% (93) 

where i/"! = '01 ('^^y, '^Cd)' J° and j'' are ac and dc current 
density values, and -R^^ are even and odd parts of the 
i?a, respectively. 

It is worth noticing that Eqs. ([89]) and (|90|) have the 
same structure as Eqs. (j46|) and (|47l) for the dc + ac re- 
sponse at e ^ 1. Actually, if we change Z^ and Z^ ^ 

in Eqs. (HH]) and ([90]) by 1/0=*= and pjj^'f , respectively, we 
obtain then Eqs. ([46|) and ([47]) . 

So all conclusions following the discussion about a 
structure of these equations can be repeated for the 
Eqs. ([89|) and ([90]). 

It is interesting also to analyze an anisotropic power 
absorption in the limit of e ^ 1 , given by Eq. ([86| in the 
previous Section IV. E. 

Let us put Gi — Gi + Gj^, where Gf are presented 
by Eqs. ([M|) and ([Mj). In case where e = 0, Gi = Gj^ = 
Gi«^,0 = where = ^cosa, 

iy = COS a and 

p = p/[l + cosa • nlmi/'f (n)/,^"], (94) 

Note, that Eq. ([M|) at a = yields p = pT(e = 0), 
where pT is given by Eq. ([88]). 



G. Linear ac response 

Here we assume that j = jdcj^jac^iut ^-^^ ^j-^g alternat- 
ing current is small (j""^ <C 1). There arc three different 
ways to derive linear (in j"'^) impedance Zu at arbitrary 
value of j'^'^. 



13 



The first way is to use general expression for 
Zi{j'^'^,j'^'^) (see Eq. (1571) ) derived by the method of ma- 
trix continued fraction at arbitrary magnitudes of the j'"^ 
and j'^'^. If we take into account that ipi{j°''' = 0) = 0, 
then it follows 



L = lim Zi = l-z(dV'i/dr=)|,.c=o. 



(95) 



This method is the most general and powerful if we can 
calculate Zi(j"^ 

The second way is to calculate Zu by means of making 
the perturbation expansion of the {r™){t) (see Eq. (ITO)) ) 
in powers of ^ 1 in the form 



(r" 



(96) 



where (r™)! = Am{i^)f''e''^* and the subscript "0" de- 
notes the statistical averages in the absence of the ac and 
the subscription " 1" the portion of the statistical aver- 
age which is linear in the ac. Whereas the {r'^)o satisfies 
Eq. the complex amplitude Am{oj) (for m > 1) can 
be presented (see details in subsection 5.5 of Ref. fi^]) in 
terms of the infinite scalar continued fraction S'm(cj) as 



Ai{oj) = 2z^(-l)" n 5™^,„(w), 



n=l 



(97) 



where 

Sm(w) 



m— 1 



1/2 



m -J n 



1/4 



ii^f I Vide I ™+l 
m+1 ^ '-^ g 



(98) 

and Sm{u} = 0) = Sm (see also Eqs. and ((Ml) ). 
Using Eq. (|97p and taking into account that A^i{uj) = 



-Al{—uj) we conclude that 



(sinx>i = (^/2)((r)l - (r-i)i) - i?(c.)j-e'' 



(99) 



where B{u}) = (i/2)[Ai(tj) + Al{-u})]. Then from ex- 
pressions for the {vx)ti, taken at j°-'^e^'^*- , Eqs. (|63|) and 
(|99p follows that dimensionless linear impedance is 



Zu{u:) = \-B{u). 



(100) 



At last, the third way to calculate the linear impedance 
gives an approximate analytical expression for Zii{uj) 
within the frames of the method of effective eigenvalue 
(see details in subsections 5.6 and 5.7 of Ref. [^'^]). Fol- 
lowing this approach we can express the dimensionless 
linear impedance in terms of the modified Bessel func- 
tions Iu{z) as 



ZiK^,g,r^) = l 



A+m(.9) 



(.9) 



1, 



2'lf,{g){\ + iLUT) I^f,{g){\* +iuJT)' 



(101) 



where 



A = /^(g)/i+^(g)/[2 / I^it)h+^{t)dt] (102) 

^0 



is an effective eigenvalue!^ and /i = igj'^'^. It follows from 
Eqs. PUT|) and pil^ that at cj = 



= l-Re[{2/lf^ig))j'l,{t)h+,{t)dt], 



(103) 



where i^o{j'^'^) is given by Eq. (|35|) . Note also that the 
right-hand side of Eq. (|103p is the exact expression for 
the dimensionless static differential resistivity in an an- 
alytical form. In the limit j'^'^ — from Eq. (|103p fol- 
lows the well-known result of Coffey and Clem^-*^ (see also 
Refs. [!^'!'^]). Actually, in this limit fi = and 



A = A*=/o(5)/i(.9)/[/o'(.9)-l]- (104) 
t/oo + (wt)^ -I- za;r(l - i/qo) 



As a result 

Zi,(c.,5,/^ = 0) = Z?, = 



1 + (ojt)^ 



(105) 

where j/qo = '^oW''^ = 0) = l//g(g) is the flux creep 
is,ctorM>lLll and 

T = f/A = f[/2(.g)-l]//o(5)/i(.g) = t{1-i.oo)[Io{9)/Ii{9)] 

(106) 

is the characteristic relaxation time. If Z^^ = p^/ ~ Kiii 
where and Cj*; are linear resistivity and reactivity in 
the absence of the dc current, respectively, then from 
Eq. (fTU5)) (see also Eq. (0)) follows that 



p?,(^,5) = l- 



1 - I'OO 
1 + (wt)2 



r° - 



UJt{1 ~ Vqq) 

' l + (wr)2 



(107) 

As expected, in the limit of zero temperature (g — >■ oo) we 
have that vqq — >■ 0, t — >■ f and the results of Gittlemann 
and Rosenblun>i2 (see also Eqs. ([T]) and ([2])) are following 
from Eqs. (fTUT)) . 



H. Nonlinear impedance and harmonics response 

Let us consider strong nonlinear effects in the ac 
impedance of a sample subjected to a pure ac drive 
dimensionless current density ^"coswt, where = 
|j°'^l/ic- In the following we will discuss the behavior 
of (^'', ri, g)-dependent impedance for simplicity in 
terms of the dimensionless ac resistivity p^^^ = pi and 
reactivity ('^'^'^ = Ci- As the angular a-anisotropy in 
these responses is omitted, the experimental observation 
of the following dependences (see Figs. 4-8) can be car- 
ried out in fact by the measurement of the and C|"t' 
responses in T geometry (see Eqs. (77), (78), and the 
definition of the Z'^j,). Figure 4 shows the dimensionless 
ac resistivity pi and reactivity (i versus ac current den- 
sity for different dimensionless frequencies $7 = ujt at 
very low temperature (g=100). 

As can be seen from the Fig. 4(a), when f2 is very 
small, the pi('C") shows several characteristic features: 
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a threshold value and a subsequent parabolic rise, 
above the threshold, with associated steplike structures. 
The threshold current density where a sudden increase in 
) starts may be defined as critical current density 
The step height decreases with increasing. The re- 
activity CilC^) shows nearly periodic dynamic 2it -jumps 
of the vortex coordinate occurring as the drive current 
density is increased (see Fig. 4(b)). The curves in 
Figs. 4(a, b) look like the similar curves discussed ear- 
lier— for the nonlinear resistance and reactance of the 
purely ac-driven resistively shunted Josephson junction 
model at T = where the overall shape and phase slips of 
these curves at several dimensionless frequencies was ex- 
plained in terms of the bifurcations in the time-dependent 
solution of the equation for the phase difference (p across 
the junctionSiiS^. Analogous bifurcations of the dimen- 
sionless coordinate x versus dimensionless time t/tt in 
our problem at T = can be calculated too. 

These bifurcations can cause sudden changes in Eyl{t) 
during one cycle of the alternating current and hence re- 
sult in steps^'^.When Q, becomes large, both the thresh- 
old and steps in pi('^") disappear and the amplitude of 
the x-jump in Ci('?") becomes larger. Also, the x-jump 
moves to large values of and the spacing in be- 
tween bifurcations becomes large which results in pi and 
Ci approaching unity. Because in our problem the abrupt 
27r-jumps of the dimensionless vortex coordinate x corre- 
spond to the overcoming by vortex of the potential bar- 



rier between two neighboring potential wells at nonzero 
temperature, our curves pi(C) a-nd Ci(^°), in compari- 
son with the curves of Refs. [^^^^^J are smoothed due to 
the influence of a flnite temperature. 

It is worth noticing that the magnitude of pi(^°|il) in 
Fig. 4(a) at < 1 is approximately equal to a constant 
which progressively increases with Q, increasing. From 
a physical viewpoint it corresponds to the enhancement 
of power absorption with the growth of O due to the 
increasing of the viscous losses accordingly to GR (see 
Introduction) mechanism. 

Now we consider the case when an ac current driven 
sample is dc current biased, i. e. the washboard pinning 
potential is tilted. In Fig. 5 we plot pi and Ci versus 
for various values of Sf" = |j'*'^|/ic at fixed dimen- 
sionless frequency f2 = 0.1 and very low temperatures 
(g=100). There are two regimes of behavior noticeable, 
corresponding to > or < S,'^- 

When £^'^ is very small (for instance, ^'^ = 0.01), the 
critical ac current density ^"(C^ — 0.01) is found to be 
equal to A ^ l(see Fig. 5(a)). When the dc drive current 
density is smaller than A (i. e. < ^'^ < A), it can be 
seen that the ac critical current density as a function 
of decreases and that both the step size and step rising 
pattern are changed. In the inset to Fig. 5(a) we plot the 
as a function of for = 0.1. A linear fit = 
-^'^ + A yields A = 1.02. Note however that this CIC'^) 
function is weakly frequency-dependent. 
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FIG. 6: The ac resistivity pi and reactivity Cj"i versus for 
various C = 0.01(1), 0.1(2), 0.5(3), 0.7(4), 1.5(5), 10(6), 5 = 
100, n = 0.1. 



When ^'^ > A, initially apparent jumps appear and 
with further increase of the ac drive current density 
the pi begins to decrease with smoothed intermittence 
steps and eventually approaches the unity. This regime 
(C^j ^'^) > is rather interesting because it shows strong 
vortex-locking effects in the ac impedance, similar to the 
Shapiro steps seen in the dc CVC's. For very large values 
of ^'^ (^"^ = 10), as expected, the effect of the microwave 
current density is negligible and the ac dimensionless re- 
sistivity pi approaches the unity. 

In the case of ac reactivity Ci, which is plotted in 
Fig. 5(b), the smoothed x-jump is not affected by the 
increase in the dc current density however, the ampli- 
tude of this jump reduces substantially for larger values 
of C''. 

In Fig. 6 we plot pi and Ci versus ^'^ for different ac 
drive current densities C° = 0.01,0.1,0.5,0.7,1.5,10. It 
can be observed from the Fig. 6(a) that for all values of 
f , as increases, pi initially increases, reaches a maxi- 
mum (very narrow and high for < 0.1) and approaches 
unity eventually. On the other hand, Ci: decreases for 
C < 0.1 slowly increases below < 1 and then sharply 
decreases, having in the vicinity of ~ 1 deep mini- 
mum, and then approaches to zero. So the occurrence of 
the x-jump can be seen clearly when ^'^ is small, whereas 
large values of either or ^'^ diminish the effect of the 
x-jumps. 

In Fig. 6(a) the pi{£,'^) dependences demonstrate sev- 
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FIG. 7: The frequency dependence of p^rj,(Q) (b) for various 
9, (I) t = 1 (solid lines), (II) = (dotted lines). 



eral main features. First, the curves, calculated at <C 1 
show the progressive shrinking of the flux creep range 
(where the pi(C'*) 1) with the increasing. If we de- 
fine the (C) the dependence of the dc critical current 
on the value of a small ac driving, then we can show that 
~ 1 — ^° at <C 1. The physical reason for such be- 
havior is obvious. Second, the appearance of a high peak 
in pi{£,'^) near = 1 for can be simply explained 

from an examination of the dc CVC curves, calculated 
at — > 0. In this limit it is evident that a dynamic dc 
resistivity (taken in the vicinity of the ^^), which equals 
to the derivative of the dc CVC with respect to the is 
strongly enhanced at T — > 0. Third, taking into account 
an analogy between Brownian motion in a tilted peri- 
odic potential and continuous phase transitions^^, one 
can say that a threshold type phase transition in the 
vortex motion along the x-axis occurs between the "lo- 
cahzed" vortex state at < S.^iS,"') and the " delocalized" 
running state at ^'^ > (^")- If consider only the lin- 
ear impedance response (i. e. Zi^l{ui) does not depend 
on ^"), this phase transition takes place at = 1 and 
at the x-point^^ where d'^Up/dx^ = 0. Forth, as it was 
shown recently in^®, a strong enhancement of the effec- 
tive diffusion coefficient D of an overdamped Brownian 
particle in a tilted washboard potential near the critical 
tilt may occur; that, in our case, D(^'^) may have a peak 
in the vicinity of = 1. 

The consequences of the D-enhancement we analyze 
with the aid of Fig. 7 where the frequency dependence 
of pi{Q\^'^ = 0) (monotonic curves) and pi{il.\£,'^ — 1) 
(nonmonotonic curves) calculated at = for three 
different temperatures {g = 10, 20, 50) are shown. The 
monotonic curves = agree with the results of Cof- 
fey and Cleni^^ who, in fact, calculated the temperature 
dependence of the depinning frequency (introduced at 
T = in^°) in a nontilted cosine pinning potential. In 
contrast to this monotonic behaviour, the nonmonotonic 
curves (C'^^l) demonstrate two characteristic features. 
First, an anomalous power absorption (pi ~ 1,6 2,8) 
at very low frequencies. Second, a deep minimum for 
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FIG. 8: The third harmonic response Re^a and \Zz\ versus 
Q. for various ^ = 0.01(1), 0.1(2), 0.3(3), 0.5(4), 0.7(5), C'' = 
0.9, = 30. 



the power absorption (pi ~ 0, 3 0, 6) at T-dependent 
^min- The appearance of this frequency- and tempera- 
ture dependent minimum at = 1 may be related to the 
resonance activated reduction of the mean escape time of 
the Brownian particle due to an oscillatory variation of 
the pinning barrier height^!. 

At last we consider the frequency dependence of 
the k-th transformation coefficient amplitude, i.e. 
\Zk\{^\j''',j°',g) taken at different values of the dc and ac 
current densities and inverse temperature. Here we point 
out only the summary of the |Z/j|(ri) curves behaviour 
because a more detailed description of these results (in- 
teresting for applications) will be published elsewhere. 

The main feature of the frequency dependence of the 
k-th harmonics is the appearance (see Fig. 8) of a pro- 
nounced maximum at ^'^ = 1 for ^Imax = 0.1 in both 
KeZk{^) and jZ^Kfi) curves. The magnitude of the 
maximum is increasing with ^"-decreasing to zero and 
(/-increasing to infinity. For example, the magnitudes of 
KeZki^lmax) and \Zk\{ftmax) are approximately equal to 
20 for fc = 3,4 at 5 = 30, = 1 and C° = O.Ol.The 
emergence of this maximum and its growth with the tem- 
perature decreasing {g — >• oo) is related to the origin of a 
singularity of the S'{u})-type in the ^-dependence of the 
linear impedance response of the overdamped shunted 
Josephson junction at T=0 and < 1^—. Note also that 
ReZfc(i7) becomes negative in the vicinity of i^max which, 



in turn, is related to the origin of a similar singularity in 
the linear Zi(j''^|ri, g) response with the Jl-increasing. 



VI. CONCLUSION 

In conclusion, the considered exactly solvable two- 
dimensional model of the vortex dynamics is of great 
interest since a very rich physics is expected from com- 
bination of a strong dc and ac driving, arbitrary value 
of the Hall effect (note, that a big Hall effect was ob- 
served in YBCO^), and the low temperature mediated 
vortex hopping (or running) in a washboard pinning po- 
tential. The obtained findings substantially generalize 
previous theoretical results in the field of the dc [2,3] 
and ac [11-13] stochastic approach to the study of the 
vortex dynamics in the washboard planar pinning poten- 
tial. Experimental realization of this model in thin-film 
geometry^ opens up a possibility for a variety of ex- 
perimental studies of directed motion of vortices under 
{dc + ac)- driving simply by measuring longitudinal and 
transverse voltages. Experimental control of a frequency 
and value of the driving forces, damping. Hall constant, 
pinning parameters and temperature can be effectively 
provided. 

While the discussion in this paper has been entirely in 
the context of nonlinear 2D pinning-mediated vortex dy- 
namics, we are aware that obtained results are generic to 
all systems with a tilted washboard potential subjected to 
an ac driving. In this sense we are conscious of that phys- 
ical explanation of our results should be supplemented by 
several new notions widely discussed. Here we mean no- 
tions of stochastic resonance^, resonance activation2£, 
noise enhanced stability^! which may be used not only 
for interpretation of our theoretical results, but on the 
contrary, the experimental verification of some predic- 
tions of these new approaches may be performed with 
the aid of the model under discussion. 

It was shown also how pronounced nonlinear effects ap- 
pear in the ac response and the linear response solutions 
are recovered from the nonlinear ac response in the weak 
ac current limit. An influence of a subcritical or over- 
critical dc current on the time-dependent stationary ac 
longitudinal and transverse resistive vortex response (on 
the frequency of an ac-driving 17) in terms of the non- 
linear impedance tensor Z and a nonlinear ac response 
at f2-harmonics are studied. New analytical formulas for 
2D temperature-dependent linear impedance tensor Zl 
in the presence of a dc current which depend on the angle 
a between the current density vector and the guiding di- 
rection of the washboard PPP are derived and analyzed. 
Influence of a-anisotropy and the Hall effect on the non- 
linear power absorption by vortices is pointed out. 

Up to now we have considered only the vortex motion 
problem. For the future experimental verification of our 
theoretical findings we should keep in mind that they 
may be applied directly only for thin-film superconduc- 
tors in the form of naturally grown (for example, in the 
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untwined a-axis oriented YBCO film22i) and artificially 
prepared washboard pinning structures^. An application 
of our results for more general cases should take into ac- 



count that they may be supplemented by consideration 
of the complex penetration length and the quasiparticle 
contribution in the way as it was made in the papersii*^. 
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